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Aims. We present a new automatic code (ARES) for determining equivalent widths of the absorption lines present in stellar spectra. We also 
' describe its use for determining fundamental spectroscopic stellar parameters. 

Methods. The code is written in C++ based on the standard method of determining EWs and is available for the community. The code 
Q ■ automates the manual procedure that the users normally carry out when using interactive routines such as the splot routine implemented in 
5-H 1 IRAK 

C/) , Results. We test the code using both simulated and real spectra with different levels of resolution and noise and comparing its measurements 
C3 • to the manual ones obtained in the standard way. The results shows a small systematic difference, always below 1.5mA. This can be explained 
' by errors in the manual measurements caused by subjective continuum determination. The code works better and faster than others tested before. 
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1. Introduction 

When analyzing the spectrum of a star, it is possible to deter- 
mine its atmospheric parameters, such as the effective tempera- 
ture (r e jf), surface gravity (logg), and the chemical abundance 
for several elements (Gonzalez & Lambert 1996; Fuhrmann 
et al. 1997). These are parameters that can be determined di- 
rectly from the spectra and can be used to determine other in- 
direct parameters. 

However, this technique can be as powerful as it is time- 
consuming. On the one hand, measuring chemical abundances 
in solar type stars often implies the measurement of accurate 
EW's for many spectral lines. On the other hand, one of the 
most common, accurate methods of determing spectroscopic 



Send offprint requests to: S. G. Sousa: sousasag@astro.up.pt 

* ARES Webpage: http://www.astro.up.pt/~sousasag/ares/ - Table 
1 is available in electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/cgi-bin/qcat ?J/A+ A/ 



stellar parameters such as r e ff, logg, and metallicity ([Fe/H]) 
is based on measuring the equivalent width (EW) of weak 
metal lines. By using them to impose excitation and ionization 
equilibrium through stellar atmosphere models, we can deter- 
mine the spectroscopic stellar parameters. Several authors (e.g. 
Gonzalez & Lambert 1996; Gonzalez et al. 2001; Sato et al. 
2003; Santos et al. 2004; da Silva et al. 2006; Randich et al. 
2006) use similar methods to present excellent results with high 
precision that are very consistent. 

The most common procedure for measuring EW is the use 
of interactive routines such as the IRAF 1 "splot" routine within 
the echelle package. In these cases it is usually necessary to 
"manually" find the line and to mark the position of the con- 
tinuum and the position of possible additional lines in the case 
of blending effects. Even the inexperienced reader can imag- 

1 IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, U.S.A. 
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ine the time needed to analyze a spectrum in this way, taking 
the number of lines for each spectrum into account to deter- 
mine the parameters with satisfactory accuracy. Multiplying 
this time by the number of stars, one can have an idea of the 
total time needed to study a large sample of stars for several 
chemical elements. 

With the growing amount of incoming stellar data (e.g. 
planet-host star sample, comparison star samples, etc.), it be- 
comes necessary to accelerate the analysis process. To over- 
come this problem, people have tried to build automated codes 
to measure EWs. In Sousa et al. (2006) we presented the re- 
sults of one promising code for this task. From several codes 
found in the literature, DAOSPEC (Pancino & Stetson in prepa- 
ration) seems to be a satisfactory code for measuring EWs au- 
tomatically. However, the user must pay particular attention to 
the choice of the input parameters for the code. Moreover, the 
code seems to underestimate the measurements by a significant 
amount. 

In this paper we present a new automatic code (ARES) that 
can measure the EWs of absorption lines. In Sect. 2 we report 
the basic idea behind the procedure and the system require- 
ments to run ARES. The ways to use the code are reported in 
Sect. 3, which refers to the input parameters and explains the 
numerical method inspired in the "hand" made measurements. 
In Sects. 4 and 5 we report the results of tests to the code both 
for synthetic spectra using the Sun as a reference and also for 
several real spectra taken with FEROS, HARPS, and UVES 
spectrographs. We also show a comparison between ARES and 
DAOSPEC EWs. In Sect. 6 we present new atomic parameters 
for a large sample of iron lines used for determining the stellar 
parameters. This sample of lines is then used to determine the 
stellar parameters and iron abundances for the Sun and several 
solar-type stars in an automatic way. The results are compared 
with the literature data. Section 7 summarizes the work pre- 
sented in this paper. 



number of lines (in case of blending effects) necessary for the 
best fit to the normalized spectrum around the line we want to 
measure. 

The procedure is to take a one-dimensional spectrum, a list 
of the spectral lines to measure and a file that contains the pa- 
rameters necessary for the computation. Then the code selects 
a region around each line, where it finds the continuum posi- 
tion. After this, it identifies the number and the position of the 
lines needed to fit the local spectrum. The fitting parameters are 
used to calculate the equivalent width of each line. The result 
is produced in a file denned by the user (see Fig. 1). 

The ARES code was written in C++, allowing any- 
one to use it without having any problems with software 
licenses. Not only is the language free, but the libraries 
that this code uses are also available without any addi- 
tional costs. The code can be obtained through the webpage 
(http://www.astro.up.pt/~sousasag/ares/), and the packages re- 
quired to compile the code are: 

- CFITSIO (http://heasarc.nasa.gov/fitsio/fitsio.html); 

- GCC - GNU Compiler Collection (http://gcc.gnu.org/); 

- GSL - GNU Scientific Library (http://www.gnu.org/soft- 
ware/gsl/); 

- plotutils (http://www.gnu.org/software/plotutils/). 

All software used to run ARES, with the exception of 
CFITSIO, can be easily installed in any linux flavor machine, 
either by using yum for Fedora Core linux or using repositories 
for Debian. The user can proceed by installing each software 
package individually by making use of the information pro- 
vided on each software web page. We want to thank the Free 
Software Foundation (FSF) for sponsoring the GNU Project, 
making it possible for anyone to have access to free functional 
software. 



2. The ARES code 



The idea behind this code is to reproduce the "manual" proce- 
dure of measuring the equivalent widths of absorption lines in 
spectra. The underlying approach is to try to "teach" the com- 
puter how to estimate the continuum position and to guess the 



ID FITS Spectra 



Line List 



\mine.opt\ 



ARES 

- Estimates local continuum 
around each line from the 
input list 

- Fit local line/lines 

- Compute EW, FWHM, ... 



lOutput Filej 
Fig. 1. Scheme for the use of ARES. 



3. Numerical method and input parameters 

A series of input parameters are needed to run ARES. These 
parameters are given in an input file and include the ID-fits file 
name (eg. IRAF format), the line-list, the output-file names, 
the wavelength domain to be used, the minimum interval be- 
tween successive lines, the wavelength interval to be consid- 
ered around each line, the minimum accepted EW of a line, cal- 
ibration parameters for either the continuum determination and 
noise control, and the display parameters. A more thorough de- 
scription of these parameters can be found in a "README" file 
that comes with the program package. Note that the input spec- 
trum should already be reduced and corrected by the Doppler 
shift. We also recommend making a previously lst-order nor- 
malization, especially for regions where the continuum has big 
changes. 

Some other more crucial parameters, needed to fit the con- 
tinuum and to smooth the spectra, are also defined here. Given 
their relevance, these are described in more detail in the next 
sections. 
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rejt = 0.99 




Fig. 2. A synthetic spectrum region with a signal-to-noise ratio around 
50. The filled line represents the fit to the local continuum through a 
2nd-order polynomial based on the points that are marked as crosses 
in the spectrum. We can see the variations induced by some different 
choises to determine the local continuum (effect of the parameter rejt, 
cf. Sect. 4.1). 



3.1. Continuum position determination 

The computation of the EW is done locally around each line 
that the user wants to measure. The parameter space is defined 
by the user in order to choose the local interval (in Angstroms) 
where the calculations will be carried out. Therefore, the con- 
tinuum position is also determined locally as is usually done 
when obtaining the EWs through the "manual" procedure. The 
continuum position is obtained iteratively by choosing points 
that are above the fit of a 2nd-order polynomial times a value 
that is defined in the parameter rejt. This important parameter 
is necessary for responding to the S/N of the spectra. In this pa- 



per we recommend some values for this parameter for different 
kinds of noisy spectra; however, since the "manual" continuum 
placing is very subjective, we decided to leave this parameter 
free for the user to choose the best value. 

Figure 2 illustrates the changes that the parameter rejt in- 
troduces in the selection of the points for the continuum deter- 
mination. This figure shows a region of a synthetic spectrum 
with an S/N around 50. Each panel is the result for a different 
value of the parameter rejt where we mark the final points cho- 
sen after all the iterations have been done for the continuum 
determination. It is easy to see that higher values for this pa- 
rameter, close to the value 1, will be advised for spectra with 
low noise. We present a case where the spectra has some signif- 
icant noise level, therefore a value between 0.94 and 0.97 will 
make a more accurate continuum determination. 

There is also an additional internal condition that tries to 
neglect possible cosmic rays/bad pixels in the spectra by taking 
out points that suffer a big change in the flux. This is done by 
checking the difference between consecutive points. If this has 
a high relative value (10%), then this point is not taken into ac- 
count for the continuum determination. In the synthetic spectra 
we introduced a bad pixel around 5501.2 A in the continuum 
region. It is easily seen that in all cases this pixel was ignored 
for the continuum determination. We stress that it is assumed 
that the spectrum used as input for the calculation should be re- 
duced as much as possible and calibrated in wavelength, with- 
out cosmic rays and with some type of normalization. 

Note that the position of the continuum, together with the fit 
of the line, can be seen if required by the user (setting plots Jlag 
= 1). Then the user has the choice of continuing to run with- 
out plots if the result is satisfactory (setting plots jlag = 0). 
Alternatively the user can stop the code and reset the parame- 
ters for a better result. 



3.2. Finding the lines to fit 

To calculate the EW it is necessary to fit the line with a func- 
tion. Absorption lines weaker than 150 mA are described very 
well by a Gaussian function. In this version of the code we only 
fit the lines with Gaussians, but this can be easily extended to 
other options if necessary in the future. 

When using interactive ("manual") routines it is necessary 
to point the approximate position of the line peak or peaks in 
the case of blending lines. To enable the computer to automat- 
ically find the peaks, we use some mathematical properties of 
the derivatives of a function. It is well known that the zeros 
of the derivative of a function give us the local minima and 
maxima. The zeros of the 2nd derivative give us the inflection 
points, i.e. the points where the function changes its concavity. 
The local maxima of the 2nd derivative will give us the cen- 
ter of the absorption lines directly. The best way to find these 
maxima is to use the zeros of the 3rd derivative of the function. 
This can be seen in the right panel of Fig. 3. It is clear that the 
maxima of the second derivative gives us a good estimation of 
the position of the lines in the original spectra even for cases of 
strong blended lines (e.g. 5500.6 A and 5500.8 A in Fig. 3). 
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Fig. 3. A simulated normalized spectral region and the respective 1st, 2nd, and 3rd derivatives. In the right panel it is easily seen that the local 
maxima of the 2nd derivative, which can be found using the zeros of the 3rd derivatives, estimates the center of the lines of the simulated 
spectrum very well. Note that using this we can see extremely blended lines that do not present zeros on the 1st derivative. The left and middle 
panels show how the parameter smoothder is useful for eliminating the noise when computating the derivatives. 



To determine the peak or peaks to fit the spectral line, it is 
necessary to obtain the first three numerical derivatives of the 
surroundings of the line profile. 

3.3. Dealing with noise 

Noise can be a problem, especially when determining numeri- 
cal derivatives where the noise propagates very fast (left panel 
of Fig. 3). To overcome this problem we make use of a nu- 
merical smoothing applied to the arrays of the derivatives with 
a boxcar average and a given width eliminating some of the 
noise. The use of such a smoothing parameter is illustrated in 
the middle panel of Fig. 3. The parameter smoothder is an inte- 
ger value corresponding to the width (in pixels) of the boxcar. 
In the code it is only applicable to the derivatives. The smooth 
result of an array of values is described by the expression 

w- 1 

Aj otherwise 

where A is the array to be smoothed, w the width of the box- 
car to be averaged, the number of array elements, and R 
the resulting smoothed array. We note that this parameter can 
be important regarding the numerical resolution of the spectra. 



From the tests we made using different spectra from FEROS 
and HARPS, which present a significant difference in spectral 
resolution and consequently also numerical resolution, the typ- 
ical value for this parameter may change between 2 and 4. 

Even using this numerical trick, noise can be a problem 
when identifying the lines. If the noise is high enough, the pro- 
cedure can mistakenly identify lines that are very close to each 
other and that is only one line in reality. To overcome this pos- 
sible problem the user can set the parameter lineresol that tells 
the computer the minimal distance between consecutive lines. 
The value is introduced in Angstroms units. A typical value for 
this parameter is 0. 1 A. 

4. Testing ARES: Using the solar spectrum 

To test ARES we performed the same tests that were presented 
in Sousa et al. (2006). We will compare the "hand" made EW 
measurements obtained using the IRAF "splot" routine to the 
ones computed by ARES. The test will be carried out using sev- 
eral solar spectra with different levels of artificial noise and res- 
olution. The original spectrum was obtained using the Kurucz 
Solar Atlas, and the noise was introduced using a Gaussian 
distribution. The artificial instrumental resolution was included 
using the "rotin3" routine in SYNSPEC 2 (Hubeny et al. 1994). 

2 http://tlusty.gsfc.nasa.gov/index.html 
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Fig. 4. ARES results for different kinds of spectra in terms of resolution and noise contamination. The larger plot shows the comparison between 
ARES EWs (y axes) and "hand made" EWs (x axes). In the box we present the slope of the linear fit to the points (dashed line), the rms, the 
number of identified lines, and the mean difference in the EW of the measurements. The filled line is a 1 : 1 line. The small plots represent the 
variations in the values presented on the box for the different spectra. On the left we fixed the resolution and changed the noise levels. On 
the right we fixed the noise level and changed the resolution. Note that these results can be improved if a more careful choice of the input 
parameters is made for each spectrum type. 



To perform the tests we chose 1 14 iron lines (blended lines 
also included) in the region [6000 A- 7000 A]. In Fig. 4 we have 
a set of plots that show the results of the tests. The large panel 
shows a comparison between the ARES EWs (y axes) and the 
"hand made measurements" (x axes) for a spectrum with S/N 
~ 100 and a resolution R = A/AA ~ 65 000 (note that R ~ 
Aclfwhm, where A c = 6500 A and/w/zm=0.1 A). In the box 
of this plot we indicate the slope of the linear fit to the points, 
the number of identified lines by ARES, the rms, and the mean 
difference of the results (in mA). As we can see from the plot, 
the points are all very close to the identity line (filled line). 
In this figure we also present several panels that illustrate how 
the values displayed in the box mentioned above change when 
fixing the spectral resolution (R ~ 65000)-left and when fixing 
the noise level (S/N ~ 100)-right. 

In these comparisons, a few points can appear to be at a sig- 
nificant distance from the identity line (2-3 points in the larger 
plot for Fig. 4). These offsets can appear for strongly blended 
lines where the fit is more difficult to compute. To check these 
cases, the user can see the log file that the code produces to 
check for bad fits in these lines. 



The slopes are always very close to 1, and the rms is al- 
ways smaller than 2 mA except in the case of low resolution. 
The number of identified lines is very high in most of the spec- 
tra and the mean difference always show a low value in the 
typical spectral type. The results are generally very good, but 
the agreement gets worst when the spectra quality decreases, 
as expected. We note that the parameters used to perform this 
test for the different kinds of spectra were chosen in a "fast" 
way, i.e. a more careful choice of the parameters for each type 
of spectra is required to obtain better results. The parameters 
are very important, and once these are adequately chosen for 
a specific type of spectrum in terms of resolution and noise 
level, the results will be equally valid under similar conditions. 
We recommend the users to make a careful selection of the pa- 
rameters, or to take the ones presented in the next section for 
reference. We also recommend the user to first use a line list 
composed with "easy" lines (isolated and well-defined) in or- 
der to see the progress of the code. The user can do this by using 
the plots provided by the code and choose the parameters that 
present the best fit to the lines and the continuum. As discussed 
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Fig. 5. Results of both ARES (right) and DAOSPEC (left) comparisons to 'hand made' measurements using the same sample of FEROS spectra. 
The dotted histogram represents the lines that are weaker than 40 mA. The thin line histogram represents the lines that are stronger than 40 mA. 



above, the continuum position is very subjective, and each user 
can choose the parameters that better fit his/her demands. 

5. Testing ARES: using real stellar data 

The code has also been applied to available spectra obtained 
with the FEROS, HARPS, and UVES spectrographs. 

5.1 . FEROS spectra; ARES vs. DAOSPEC 

A comparison to the results obtained by DAOSPEC for the 
same FEROS sample of spectra is also included. The spectra 
were obtain to characterize the planet-host stars and compari- 
son samples of stars that was not found any planet. For details 
on these spectra, see the work of Santos et al. (2004, 2005) and 
Sousa et al. (2006). 

In Fig. 5 we show the results of both DAOSPEC (left) and 
ARES (right) versus the "manual" measurements. We also plot 
a histogram representative of the difference between the mea- 
surements to give a better idea of the dispersion and mean value 
of the automatic measurements. The same FEROS spectra and 
the same line list presented in Sousa et al. (2006) have been 
used. We present here the direct result from DAOSPEC without 
any calibration, details for the computation using DAOSPEC 
can be found in Sousa et al. (2006). We also present the rms 



and the mean difference of the measurements to give a better 
idea of the result of the comparison. 

In the case of the ARES code, since the sample is composed 
of high-quality spectra, with S/N > 200 and R ~ 50000, the fol- 
lowing main parameters were used: smoothder-A, space-2 A, 
rejt =0.996, lineresol-Q. 1 A, and miniline-2mk. Using this pa- 
rameters we derived EWs that show a small mean difference of 
about 1.5 mA with respect to the "hand" made measurements. 

We suspect that part of the difference mentioned above is 
due to the fact that the spectra have different noise levels and 
that we are using the same parameters to determine EWs for all 
of them. This will also contribute to the dispersion presented in 
the plots. We cannot exclude that the subjective manual deter- 
mination of the EWs may also have a significant weight in the 
difference presented here. 

It can also be seen in Fig. 5, as found before, that 
DAOSPEC underestimates EWs for most of the measurements. 
In this aspect, ARES reacts better than DAOSPEC and presents 
a sharper distribution, with the mean difference closer to zero. 
Note that the weaker lines are the ones that present higher rela- 
tive error and are the ones that contribute the most to the small 
bump observed to the left in the histograms. This can be seen 
in the histogram that represents lines weaker than 40 mA. 

Another point that favors ARES over DAOSPEC is the 
computation velocity. While DAOSPEC took a few hours to 
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Fig. 6. Results of ARES for a sample of HARPS spectra comparing Fig. 7. Results of ARES for a sample of UVES spectra comparing the 



the value of the EW obtained by ARES (EW ARES ) with the "manual" 
value from IRAF (EW IRA f). 



value of the EW obtained by ARES (EW ARES ) with the "manual" value 
from IRAF (EW, RAF ). 



obtain the measurements, since we had to do it taking sev- 
eral small intervals for each spectra, ARES took the mea- 
surements in two minutes using the same machine. Moreover, 
ARES was able to identify more lines in the spectra sample 
fhanDAOSPEC. 

5.2. HARPS and UVES spectra 

In Fig. 6 we show the result for a sample of HARPS spectra. 
The main parameters used here were: smoothder-A, space=2 
A, rejt-0.993, lineresol-Q. 1 A, and miniline-2mA. The results 
are very promising for all of the spectra, despite the significant 
variation in the noise level of the spectra in this sample (S/N ~ 
150-250 and R~ 115000). 

In Fig. 7 we show the result for a small sample of 
four UVES spectra. The main parameters used here were: 
smoothder-A, space-2 A, rejt =0.999, lineresol-Q.\ A and 
miniline=2mA. The four spectra, used here to test the code on 
UVES spectra, were of excellent quality with a very high S/N > 
200 and R ~ 100000. Therefore the value of the rejt parameter 
is also high. 

The code seems to work very nicely with these spectra, us- 
ing the parameters presented here. 



6. Stellar parameters 

6.1. Using ARES EWs 

We re-determined the stellar parameters for the FEROS sample 
of stars presented in Sousa et al. (2006) to compare the results. 
The EWs were measured with ARES, with the same parame- 
ters as described before for this FEROS sample. After this, the 
same procedure was used, using the same iron line list for a 
spectroscopic analysis done in LTE using the 2002 version of 
the code MOOG 3 (Sneden 1973) and a grid of Kurucz Atlas 
plane-parallel model atmospheres (Kurucz 1993). 

Figure 8 shows the comparison of the parameters obtained 
using ARES versus the parameters obtained using the manual 
EW measurements. We can see from the figure that the pa- 
rameters are all very close to the identity line. The effective 
temperature has a mean difference of 22 K and the metallic - 
ity has a mean difference of 0.03 dex. These differences can 
be explained by the small difference shown before of 1.5 mA 
presented in the EWs measurements (see Sect. 6.1). 



3 The source code of MOOG2002 can be downloaded at 
http : //verdi . as .utexas .edu/moog .html 
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Fig. 8. Comparison of the manual obtained spectroscopic parameters to the obtained automatically through ARES code 



Table 1. Sample Table: Atomic parameters and measured solar equiv- 
alent widths for Fe n and Fe i lines. 
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-2.638 


41.3 


4139.93 


0.99 


-3.423 


80.4 



6.2. Extended list of iron lines 

We built a code that presents promising results in the automatic 
EW determination. This was the first step in the quest for a total 
automatic procedure to determine spectroscopic stellar param- 
eters. The second step is to determine a way to obtain the stellar 
parameters with high accuracy. The idea is to use a significant 
number of iron lines to increase the statistical strength of the 
determined stellar parameters. 

We follow the procedure presented in Santos et al. (2004) 
for derivating stellar parameters and metallicities. It is based 
on the EW of 39 Fei and 16 Fen weak lines, by imposing 



excitation and ionization equilibrium. Presently the number of 
lines are relatively small since users are restricted to a few well- 
defined iron lines in the stellar spectrum. The idea is to increase 
this number of lines, giving the possibility of decreasing the er- 
rors involved in the determination of the stellar parameters. 

Using the Vienna Atomic Line Data-base (VALD: 
Piskunov et al. (1995); Kupka et al. (1999); Ryabchikova et al. 
(1999) - http://www.astro.uu.se/htbin/vald), we requested all 
the Fe i and Fe n from 4000A to 9000A. Then using the Kurucz 
solar atlas, we searched all the iron lines allowing a good mea- 
surement. Then the log gf values were obtained. For this a 
differential analysis to the Sun was used; i.e. the log gf val- 
ues were determined using the lines present in the Sun. We 
computed from an inverted solar analysis using EW measured 
from the Kurucz Solar Atlas and a Kurucz grid model for the 
Sun (Kurucz 1993) having T eff = 5777/T, logg = 4.44 dex, 
% t = 1.00 kms~ l , and log(Fe) = 7.47. The new values of log gf 
are presented in Table 1 . 

As a test, we computed the solar parameters and 
iron abundances based on iron EW measured using a 
solar spectrum taken with the HARPS spectrograph us- 
ing the Ceres asteroid (Collection of HARPS solar spec- 
tra - http://www.ls.eso.org/lasilla/sciops/3p6/harps/monitoring/ 
sun. html). The resulting parameters were T e ff = 5762 + 12 , 
logg = 4.44 + 0.08, & = 0.94 + 0.01, and [Fe/H] = 0.0 + 0.02, 
close to the "expected" solution. 
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We also show results for a few stars in the FEROS and 
HARPS sample. The stellar parameters obtained for these stars 
were obtained using ARES to measure EW with this new line 
list. The procedure is the same as for the short line list and the 
stellar parameters are presented in Table 2 where we can com- 
pare them to previous results added from the literature. The re- 
sult of imposing the excitation and ionization equilibrium can 
be seen in Fig. 9 for the HD20201 star. The results are compati- 
ble with the previous ones obtained using the short line list, but 
since we have more lines, this is statistically more significant. 



HD2020': 




3.0 [ 




-6.0 



-5.5 -5.0 
log (RW) 



-4.5 



Fig. 9. Results for the fit of the abundances for determining the stellar 
parameters for the FEROS spectra of HD20201. The crosses represent 
the Fe i, and the filled triangles the Fe n. 



6.3. A more accurate analysis 

An alternative procedure that can improve the accuracy of the 
stellar parameters is being considered. A differential spectro- 
scopic analysis (Laws & Gonzalez 2001) between the Sun and 
a solar-type star can eliminate intrinsic errors in the determi- 
nation of the stellar parameters. If we use the same spectro- 
graph to obtain both the Sun and the target star spectra, and 
also use the same method to determine the EWs, it is possible 
that, when comparing the individual line abundances for the 
both star, intrinsic errors either on the EW measurements or in 
the log gf values can be eliminated. Using ARES for the EW 
measurements, we are now able to apply this procedure to a 
large sample of lines. 

Another alternative option that can also improve the analy- 
sis is to compute new log gf values based on ARES measure- 
ments for the solar spectrum obtained with the same spectro- 
graph (same instrumental and numerical resolution) and noise 
levels. In this way, a systematic error in the EW measurement 
made by ARES will be taken into acount in the respective new 
solar log gf value. 



7. Summary 

We present ARES, a new automated code for measuring EWs. 
We describe the method used by this new code and also the 
requirements for its use. 

The results of the test performed with ARES concludes 
that this code seems to perform very well for several types of 
spectra. There is a small difference between ARES and man- 
ual measurements that can in part be explained by the intrisic 
errors on the subjective determination of the continuum level. 
The general results of the EW measurements using real data 
show excellent results, especially in high-resolution spectra. 

We also present a new large line list of iron lines for deter- 
mining the stellar parameters to be used in an automatic proce- 
dure. The result for the Sun and for the stars of the FEROS and 
HARPS samples shows that the combination of ARES with this 
large line list can be a powerful tool for the statistical precision 
of determining automated stellar parameter. 
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Table 2. Stellar parameters for some test stars 



HD number 


T eff [K] 


logg spec [cms 2 ] 


6 [kms- 1 ] 


[Fe/H] 


N(Fei,Fen) 


<r(Fe i,Fe n) 


Ref.t 


FEROS spectra: 


HD 20201 


6064 ± 34 


4.43 ± 0.02 


1.17 ±0.05 


0.19 ±0.05 


38,15 


0.04, 0.06 


[1] 




6084 ± 18 


4.57 ± 0.03 


1.16 ± 0.02 


0.16 ±0.01 


350,57 


0.05, 0.09 


[2] 


HD 33214 


5180 ± 74 


4.40 ±0.11 


1.10 ± 0.10 


0.17 ±0.08 


39,14 


0.08, 0.13 


[1] 




5265 ± 31 


4.48 ± 0.08 


1.10 ±0.06 


0.19 ±0.02 


348,51 


0.09, 0.24 


[2] 


HD 36553 


6103 ± 46 


3.85 ± 0.03 


1.61 ± 0.07 


0.41 ± 0.06 


37,15 


0.05, 0.09 


[1] 




6153 ± 23 


4.06 ± 0.04 


1.67 ± 0.02 


0.39 ± 0.02 


349,55 


0.07, 0.12 


[2] 


HD7570 


6198 ± 39 


3.49 ± 0.02 


1.40 ± 0.07 


0.24 ± 0.05 


38,15 


0.04, 0.06 


[1] 




6204 ± 16 


4.55 ± 0.03 


1.31 ± 0.02 


0.21 ±0.01 


350,59 


0.05, 0.10 


[2] 


HD7727 


6131 ± 41 


3.34 ± 0.02 


1.18 ±0.07 


0.16 ±0.05 


38,15 


0.05, 0.07 


[1] 




6169 ± 18 


4.55 ± 0.03 


1.26 ± 0.02 


0.12 ±0.01 


327,57 


0.05,0.10 


[2] 


HARPS spectra: 


HD 1021 17 


5672 ± 22 


4.27 ± 0.03 


1.05 ± 0.02 


0.30 ± 0.02 


37,8 


0.03, 0.03 


[3] 




5671 ± 18 


4.24 ± 0.03 


1.06 ± 0.02 


0.28 ± 0.01 


286,56 


0.05, 0.10 


[2] 


HD 212301 


6256 ± 28 


4.52 ± 0.02 


1.43 ± 0.06 


0.18 ±0.02 


36,9 


0.03, 0.01 


[1] 




6273 ± 21 


4.64 ± 0.04 


1.32 ± 0.02 


0.16 ±0.02 


278,53 


0.05, 0.13 


[2] 


HD4308 


5685 ± 12 


4.49 ± 0.02 


1.08 ± 0.03 


-0.31 ±0.01 


35, 8 


0.01, 0.02 


[1] 




5675 ± 13 


4.43 ± 0.02 


0.91 ± 0.02 


-0.33 ± 0.02 


293,56 


0.04, 0.08 


[2] 


HD 69830 


5385 ± 19 


4.37 ± 0.03 


0.80 ± 0.03 


-0.05 ± 0.01 


35, 8 


0.03, 0.02 


[1] 




5410 ± 17 


4.35 ± 0.04 


0.83 ± 0.02 


-0.07 ± 0.01 


284,55 


0.05, 0.15 


[2] 



t The instruments used to obtain the spectra were: [1] Sousa et al. (2006), [2] this paper, [3] Santos et al. (2005); 
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